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Abstract. - An outstanding question in the physics of soft jammed packings concerns the nature 
of the correlations that arise near the unjamming transition. In this work, we treat unjamming 
as a constraint satisfaction problem and demonstrate that a static correlation function, which 
probes sensitivity to boundary conditions, exhibits a diverging correlation length as the packing 
is decompressed. This length scale is related to isostaticity, and has been connected earlier to the 
existence of soft modes in isostatic packings. The form of the correlation function is remarkably 
similar to one predicted in the mosaic theory of the glass transition, and we relate the nature of 
the growing correlations to the entropy of packings. 
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Introduction. — There is a class of systems in soft- 
dondensed matter, which exhibit pronounced anomalies 
in their dynamics without any dramatic changes in their 
static properties. The viscosity of supercooled liquids in- 
creases dramatically as the temperature is lowered result- 
ihg in the liquid falling out of equilibrium at a nominal 
".glass transition" temperature [S]. Systems that interact 
V^ia purely repulsive, short-range forces, such as granular 
materials can jam or unjam in response to external driv- 
ihg. The dynamics in these systems become increasingly 
heterogeneous as the non-flowing state is approached, and 
there are remarkable similarities between the dynamical 
correlation functions of these two systems [Tallin] , which 
clre strikingly different in any standard measure. Granular 
materials are athermal, and have frictional forces whereas 
supercooled liquids are in thermal equilibrium until they 
fall out of equilibrium at the glass transition temperature. 
The one feature that connects these systems is the exis- 
tence of many amorphous mctastable states [T7] . 

In mean-field theories of the glass transition such as 
the random-first-order theory, a non- vanishing complex- 
ity or configurational entropy does lead to the definition 
of a static, point-to-set (PTS) correlation length that 
probes the growing influence of the boundary [T^[T^ . 
The divergence of a PTS correlation length would indi- 
cate a true thermodynamic phase transition, where by 



deflnition boundary conditions affect the interior com- 
pletely. A question that naturally arises is whether such 
a PTS measure can be used to detect phase transitions 
in athermal systems, and in particular whether the jam- 
ming/unjamming transition in granular systems can be 
characterized by a similar correlation length. In addition 
to its application to the glass transition problem, the PTS 
correlation function has been developed and widely used to 
monitor the structure of the solution space of constraint 
satisfaction problems involving discrete or boolean vari- 
ables such as those found in coloring problems and satis- 
fiability problems, for instance, XORSAT [H]. 

Jamming and unjamming can both be viewed as "con- 
straint satisfaction" problems. If grains are idealized as 
being infinitely rigid, then the constraints for jamming are 
zero overlap between grains, and the degrees of freedom 
are the positions of the grains. There has been consid- 
erable recent work analyzing the solution space of hard 
sphere packings as one approaches jamming, and there 
is evidence for a diverging length scale associated with 
this entropy- vanishing transition [M]. In systems inter- 
acting via soft repulsive potentials cohesion arises only as 
a consequence of compression, and a solid can unjam (fail 
mechanically) with all elastic moduli vanishing as the as- 
sembly of grains is decompressed [TO]. For this type of 
unjamming the constraints are the equations of mechani- 
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cal equilibrium, and the degrees of freedom are the posi- 
tions of the grains, and the contact forces. Experiments in 
granular systems have shown that unjamming through de- 
compression is accompanied by large fluctuations of stress 
PiriS]. It has been argued that for frictionless grains, the 
jammed configurations at this unjamming point, referred 
to as Point J |19| are critical with diverging correlation 
lengths, and soft modes PUZ]. So far, however, it has not 
been possible to measure a diverging correlation length 
from measurements of any static correlation functions |10| . 
The existence of a diverging length scale has been inferred 
from characteristics of the soft modes [TU] . In this work, 
we demonstrate that the unjamming transition in two- 
dimensional assemblies of frictionless disks is character- 
ized by a diverging length in a PTS correlation function, 
and that this correlation function is intimately related to 
the entropy of jammed, mechanically stable states. 

Model. — In contrast to the constraint satisfaction 
problems that have been analyzed using the PTS correla- 
tion function, jamming/unjamming poses new challenges. 
The degrees of freedom involved in unjamming are con- 
tact forces and positions of grains, which are continuous 
degrees of freedom. Given the complexity of the prob- 
lem, we study a minimal model of static granular packings 
known as the force- network-ensemble (FNE) [T] , which ex- 
plicitly incorporates local force balance. For hard but not 
infinitely rigid grains, small deformations can lead to large 
changes in the force magnitudes. The length scales govern- 
ing packing geometries and contact forces are, therefore, 
well separated [T]. In this limit, the jammed system can 
be modeled by the FNE, which is a flat distribution of 
all forces that satisfy the constraints of mechanical equi- 
librium on a given geometry [T]. The unjamming prob- 
lem then reduces to a constraint satisfaction problem on a 
particular random network, where the degrees of freedom 
are the force magnitudes. The FNE is a minimal model, 
which has been shown to capture the robust features of 
contact force distributions, and mechanical properties of 
static granular packings. Since we are interested in explor- 
ing correlations that are strong indicators of unjamming, 
and not sensitive to specific force laws, we study the PTS 
correlation function in FNE. We measure PTS correlations 
as a function pressure (allowing the set of quenched ge- 
ometries to change with changing pressure) and find that 
a static length scale does diverge as the system unjams. 
The origins of this length scale can be understood from a 
simple bulk-surface argument. In addition, we explore the 
structure of the solution space of the constraint satisfac- 
tion problem, and discuss its implications for the entropy 
of granular packings. Parenthetically, the concept of the 
force-network ensemble is on firmer footing for frictional 
rather than frictionless grains since the indeterminacy ex- 
ists even for infinitely rigid grains |20] . and some of our 
conclusions could carry over to frictional systems. 

To study unjamming, we need to sample a large set of 
amorphous geometries, and construct the corresponding 



FNEs. We use the O'Hern protocol [3] to generate an un- 
biased sampling of just-touching disk packings. Packings 
at arbitrary overcompressions are then generated by in- 
creasing the packing fraction in a controlled manner. In 
this work, we have only studied packings obtained from 
a harmonic interaction. The geometries obtained via this 
protocol are then used to construct the FNE. For a given 
geometry, PTS is well defined as the overlap between force 
networks, and can then be averaged over various realiza- 
tions of the quenched geometry. 

In a d— dimensional granular assembly of M grains and 
z contacts there are | contact force variables, which are 
constrained by dM mechanical equilibrium (ME) equa- 
tions. A minimum condition for finding valid (i.e. satis- 
fying ME) force networks is that | > dM. For overcom- 
pressed packings that satisfy this condition, FNE assumes 
a flat measure for the set of force networks /, which satisfy 
the matrix equation 

Af = b (1) 

It is the nature of this solution space that will determine 
the structure of the PTS correlation function. Through- 
out this work, we will refer to the vector / as a "force 
network," so as to avoid confusion with "force vector," 
which one typically associates with a particular force act- 
ing on a body. The force network is a list of contact force 
magnitudes. Eq. ((TJ expresses the ME constraints for 
each grain in the packing. / is a vector containing z con- 
tact force magnitudes and ^ is a matrix which contains 
all of the information about the geometry, and extra rows 
that generate any global constraints on the force network 
[T] . The values of these global constraints are represented 
by b. In our work, we choose to flx all three independent 
elements of the stress tensor CTq./? = J2i=i J2j>i ^ijfij' ^° 
that (for 2D packings) there are three global constraints 
expressed in A and three non-zero constants appear in 6. 
Constructed in this way, A has z columns and 2M+3 rows. 
When A is rectangular with a number columns larger than 
the number of rows, it will have a null space of dimension 
z — 2M + 3, if the remaining 2M + 3 rows arc all indepen- 
dent (this would not be true for crystalline packings, but 
will generically be true for amorphous packings). 

Sampling Force Network Solutions. The null 

space basis vectors {g} span a space of solutions for the 
homogeneous equation Af = 0, so that these solutions are 
at (Tq/j = 0. a real packing has non-zero ctq^, so one needs 
a particular solution /o at the fixed ctq,^ . All /o -I- eg with 
an arbitrary coefficient c are valid force networks. We 
obtain /o from the force law used to form the quenched 
geometry. Finally, while /o + eg is a solution to eq. [TJ it is 
not generally a solution that applies to a granular packing 
because the forces are not generally positive, and the posi- 
tivity constraint has to be enforced separately. We sample 
solutions by constructing a random walk in the null space 
of A. Assuming each valid solution is equally likely, a new 
solution /' is found by choosing a random step size c from 
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a uniform distribution, as well as a random direction in 
the null space g and adding eg to a previous solution /: 
/' = / + eg. In this way the sampling method is a random 
walk in the null space starting at /q. The random walk 
is subject to reflecting boundary conditions resulting from 
the positivity constraint on the individual forces as well 
as the global constraints built into A, including the ones 
coming from the boundary that are needed to define the 
PTS correlation function. In earlier studies of FNE [2] 
a "wheel move" has been used to sample valid force net- 
works. In contrast to the wheel move, the moves presented 
here are intrinsically non-local. The null space basis vec- 
tors have of order z non-zero components, so that each 
move f ^ f + eg will generally change all, or most, of the 
contact forces simultaneously. These non-local moves obey 
detailed balance, and are subject only to the constraints 
that arc explicitly built into A. We have checked that 
our sampling method reproduces single force distributions 
for real geometrically distinct packings by comparing the 
single force distributions of such packings to those of the 



FNE(see figure 1(b)), especially close to the unjamming 
transition. 

Point-To-Set as a Measure of Correlations. 

The PTS correlation function provides a measure of how 
far boundary effects penetrate into the interior of a sys- 
tem [Slinill], and is thus an indicator of a critical point. 
Generically, one creates a reference configuration of ran- 
dom variables, such as the spins in the p-spin model, or 
the positions of particles in an inherent structure of the 
Lennard- Jones system that satisfy the constraints. A new 
configuration, satisfying the constraints, is then created 
inside of some boundary B{r) while outside of the bound- 
ary the configuration is kept "frozen" , and an overlap func- 
tion is defined to measure the covariance between the old 
and new configurations. In the Lennard- Jones (LJ) sys- 
tem, for example, the PTS correlation function q{R) j6| is 
defined as the overlap between density functions of the ref- 
erence state and the "pinned state," generated by freezing 
the L J particles outside of a boundary and allowing the dy- 
namics to continue inside of the boundary. The PTS corre- 
lation function in the LJ system exhibits a growing length 
scale, indicating an increasing susceptibility to boundary 
conditions as the liquid is cooled towards the glass tran- 
sition 0. The jammed (mechanically stable) states of 
frictionless disks are inherent structures of a system with 
short-range repulsive interactions obtained by quenching 
from infinite temperature. Keeping in mind that the only 
source of cohesion in these materials is imposed pressure, 
we would like to study the behavior of the PTS as the 
pressure on the packings is reduced. To construct the 
PTS correlation function for the unjamming transition in 
frictionless disks, we therefore study packings at a given 
pressure, where a region outside of a given boundary is 
held fixed while the interior region is allowed to change in 
a way that samples new granular packings without chang- 
ing the pressure. If unjamming through decompression is 



indeed controlled by a critical point, then the PTS cor- 
relation function should exhibit a growing length scale as 
the pressure is reduced to zero. 

There is an inherent difficulty in constructing overlap 
functions for granular packings with different geometries, 
and in principle different topologies. Within the FNE 
model, we avoid this difficulty by fixing the geometry 
and only allowing the contact forces to change to calcu- 
late the PTS correlation function for a given geometry. 
The geometry-specific PTS correlation function {C{R)) is 
computed by building f{R) from /o with all of the con- 
tact forces outside of a radius R in /o left unchanged, 
but the contact forces inside of R allowed to change in 
a way that satisfies mechanical equilibrium. The overlap, 
C{R) = fo- f{R), where /o includes only grains in the core 
of the region, away from the boundary. Since they share 
a single underlying geometry, the two force networks are 
of the same length. The magnitudes of the core of each 
network are divided out so that the largest value of C is 1. 
ME equations associated with grains which share contacts 
that are inside of R are included in the set that define A, 
however, these particular equations cannot be set to zero. 
Instead, these equations are set equal to the non-zero net 
forces which are the sum of all the frozen forces associated 
with the boundary grains. Two force networks are never 
orthogonal to each other because of the positivity of the 
contact forces. But, if we ignore the positivity constraint, 
the overlap of two completely uncorrelated unit vectors 
should simply be the projection of a randomly chosen unit 
vector onto one of the axes in a z-dimensional space. There 
are z components to a randomly chosen unit vector, and 
each component squared is roughly proportional to 1/z, 
so that the overlap of two completely uncorrelated force 

networks would be oc 1/w -^z, where Mc is the number of 
grains within the core. We expect that C{R) will approach 
this value asymptotically at large R. The overlap C{R) 
is averaged over all f{R) sampled by the random walk in 
the solution space with the imposed boundary constraints 
to obtain {C{R)). To obtain the PTS correlation function 
as a function of pressure, {C{R)) is then averaged over 
many geometries at a given pressure [3] . Throughout this 
work, we will use () to denote an average with respect to 
the FNE for a fixed geometry, while Qg will denote an 
average over many geometries at a given overcompression. 

Length Scales and Scaling. — For a series of sys- 
tem sizes ranging from M = 30 to AI ~ 900 grains, and 



overcompressions ranging from 



10"-^ to 10" 



measurcd {C{R)). For each geometry, there is a unique 
scale, po = Ro/d, below which {C{p)) ~ 1, precisely, since 
the bounding region is small enough that the constraints 
on the contacts which are interior to the region cannot 
be satisfied except by the original network, /o, and there 
are no other solutions to the ME equations. For R > Rg 
{p > po), the correlation function decays monotonically 
to its asymptotic value, C'{L/d). Singular value decompo- 
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(a) 



(b) 



Fig. 1: (a) The measurement of C{R) "freezing" the contact 
forces on the boundary of a box of size R and then sampling, 
(b) (sohd lines) The single contact force distributions assem- 
bled from 40 different 900 grain geometries at (red) 5(j> = 0.1, 
(green) 5(j) = 0.05, and 5(f> = 0.01. (dashed lines) The single 
contact force distributions of the FNE sampled using a single 
amorphous geometry at (red) 5(f) = 0.1 and (blue) 5(f) = 0.01. 



sition of the matrix A can be used to identify po ^-s the 
scaled bounding ball size at which the measured nullity 
becomes greater than 0. Fig. [2] shows poj and its value 
averaged over all of the geometries at a given pressure and 
for a range of system sizes (measured in terms of number of 
grains M). For the largest system size, a fit of (po)^ versus 
the {P)g yields an exponent of v =0.461±0.012. The inset 
of fig. [21 shows finite-size scaling results for (/9o)g- We ob- 
tained data collapse over five decades using the finite-size 
scaling form -^S^ = g[L^/'^{P)g) based on po diverging as 
P~^, and taking L oc M^/^. We observe the best collapse 
for V = 0.46, with a visible failure of the collapse outside 
of the interval ±0.02. The figure shows only the region 
close to L^^'^{P)g ~ 1, The results shown in Fig. [5] pro- 
vide compelling evidence of a diverging length scale that is 
related to "ordering" in force-network space. The data for 
{{C {p)))g for different system sizes show that for p > pQ, 
increasing the system size simply continues {{C{p)))g to 
lower values with no other length scale setting in. Figl3] 
shows (C {p))g, measured as a function of p = -f, the di- 
mensionless length scale in terms of the grain diameter d. 
As seen from the figure, for p = -^ << {pojg, a pressure 
dependent scale, the correlation function is identically 1. 
The form of {{C{p)))g reflects the characteristics of the 
relevant force network solution space, and can be modeled 
as [S]: 

where Pr,f (p) is the probability of finding a force network 

different from fo, Pf^{p) = 1 -Pf^f^ip), and go = (/ ■ /o) 
is the average overlap of two unequal unit vectors that 
are constrained by positivity, and represents the "uncor- 
related" limit of C{p). A connected correlation function 
can, therefore, be defined as C{p) — go = Pr (p)(l ~ ^o)- 
Because the number of components of the unit vectors 
change with pressure, qo is a function of pressure. The 
connected correlation function obtained by approximat- 
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Fig. 2: The critical subregion size po is plotted against pres- 
sure for each of our geometries with 900 grains. (Inset) Finite 
size scaling is used to collapse po over the range of pressures 
for 150,300,450,600,750, and 900 grains, using the exponent 
u = 0.461. Different colors (symbols) correspond to different 
overcompressions, with the open symbols ranging from 0.01 to 
0.1 and the filled symbols ranging from 0.001 to 0.008. This 
symbol legend corresponds to the other figures of this article 
as well. 



ing qo by C{p = L/d) shows a rapid decay to zero and a 
collapse for different pressures after scaling p by po (Fig. 
[31). 

To compare the predictions of the model, Eq. [51 to the 
measured C{p), we need to estimate both pr^r (p), and 
go- As argued earlier, the overlap of two uncorrelated unit 
force networks depends on the number of components of 
the force network. Also, in the isostatic limit z — > zoj 
go ^- 1 and therefore we expect go — A/^f^, where z is 
the total number of contacts for a given pressure and using 
Mc oc (f). The probability of finding a force network differ- 
ent from /o depends on the volume of the solution space 
(y (p)) of ME force networks for a given bounding ball size. 
The "force network space" is a high dimensional space (of 
dimension equal to the number of contacts) and is sparsely 
populated by valid force network solutions. Embedded in 
the force network space is a solution space which is the 
null space of A shifted to an origin defined by/o, and has 
a dimension equal to the number of contacts in excess of 
the ones needed to satisfy the constraints of ME. This 
space is dense with solutions, since it is spanned by the 
null vectors of A. In addition, since each null vector is an 
equally likely component of a solution, the space is roughly 
(hyper)-spherical. Since the space is dense and spherical, 
the number of solutions for a given geometry and pres- 
sure goes as the volume of a hypersphere with dimension 
equal to the number of excess contacts, Sn. The probabil- 
ity of finding a force network different from /o is, there- 
fore, given by Pf^f^{p) = (V(5n(p)) - l)lV{5n{p)), where 
Sn{p) is the average number of excess contacts for a region 
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which is nonnegative for p > po ~ 
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Fig. 3: (Inset) The measured PTS correlation function, with 
different symbols corresponding to different values of < P >g. 
The highest pressures correspond to the curves in the inset 
which decay to the lowest values. The dashed lines are the 
predictions from eq. [2] using the measured values of 5n (p). In 
the main figure, the connected correlation function is plotted, 
as well as the model prediction for the connected part. 



of size p, and we have explicitly subtracted the one coming 
from the single solution that is identical to /g. The con- 
nected correlation function should, therefore, be given by: 
C(/3) — 9o = v(Snp) ■ ^^ ^S- ISl'^6 compare C{p) obtained 
from the model (cq. [2]), using measured values of 6n{p), 
and V{Sn{p) (discussed below) to the measured values of 
C{p). The figure also compares the connected correlation 
functions, which has a shape that is very similar to that 
predicted by the Random First Order Transition (RFOT) 
theory of glasses [T31I1I] . 

Solution Space Volume and a Mean Field Model. 

— In order to understand the origin of the {po)g behav- 
ior with pressure, we construct a meanfield model based 
on the isocounting argument discussed earlier, but now 
including a boundary contribution [4j to the number of 
constraints. If Sn (R) is the number extra contacts within 
a ball of size R beyond the number of constraints (or the 
nullity of the geometry matrix), then Sn (R) — ^^^ — B, 
where M is the number of grains inside of the ball, 
Sz = z — Ziso averaged over the packing, and B is some 
boundary contribution. Assuming homogeneity, the num- 
ber of grains inside the ball is related to the overall pack- 



ing fraction: 



A/ 71 



zw 



M 



cR' 



Ac(f)p'^, where 



c = 72 — 1 ,\ is the bidispersity factor, which has been 

kept fixed throughout this study, with a = 1.4. The 
boundary term can be estimated by first assuming that 
along the border of the bounding box, there are ^/M 
grains. Since the box has 4 sides, and there are 2 con- 
straints per grain: i? = 4 • 2 • a/M = 16-y/c0p, so that the 
meanfield expression for the nullity is: 



size po for which C [p) first becomes less than 1 is. there- 
fore, —jS= — . A length scaling of the form 4- has been ob- 

tained earlier through an analysis of soft modes. [3. Since 
it is known that, on average, the deviation from isostatic- 
ity scales with the pressure p as 5z (x. p^'^, the meanfield 
argument would predict that po diverges as p — > with an 
exponent of 1/2. This length-scale exponent has been de- 
duced earlier from a crossover behavior of the vibrational 
density of states in simulations of soft- repulsive grains |22| , 
and the meanfield argument presented here is similar to 
the argument used to deduce a diverging length scale as- 
sociated with isostaticity [7] . The numerical value of the 
exponent deduced from the PTS correlation function is 
close to the meanfield value. 

By its very definition, the PTS correlation function 
probes the size of the space of solutions that satisfy a 
given set of constraints. In light of the similarity of our 
algorithm to a random walk in a high dimensional space, 
we calculated a radius of gyration, defined as the "dis- 
tance" in the solution space from /o to the new solution: 

i?^ = ( / — /o 1 .If the solution space is roughly spherical, 

this gives a measure of the linear size of the solution space 
for a given quenched geometry and bounding ball size. If 
the walk is in a high dimensional space, it will spend most 
of its time exploring the surface of the solution space. Af- 
ter some initial equilibration time required for the random 
walker to reach the surface of the solution space, there 
is some solution fmax which is "farthest away" from the 
origin /o, and if the solution space is roughly spherical, 
I /max I is left unchanged. Taking the asymptotic value of 

\jmax\j J^g ^ \Jmax /o I ^ [Z^i=l^i9ij ^ Z^i=l '^i i 

which can be approximated as Sn[f]'^, where [/] is the av- 
erage single contact force and the square brackets are used 
for spatial averages over a particular geometry. From this 
result we see that i?g/[/]^ can provide a measurement of 
Sn (p) . The functional form of Sn (p) suggests a scaling 

form that is given hy Sn ( — ) = j^G [ — ), with G having 

a quadratic form. This scaling is found to collapse the 
Rl data well, assuming Sn{p) = i?^/[/]2 (fig. g]). The 
master curve found from this scaling indeed is a quadratic 
function of p. In (fig. 21), The mean field equivalence be- 
tween Sz and pressure is used. The collapse illustrates 
the robustness of the scaling form for Sn as well as the 
approximation of R^ as Sn[f]^. 

The model for C{p) (Eq. ^ involves V {Sn{p)), the 
volume of the solution space of force networks for a given 
p. Assuming the space to be a hypersphere of radius Rg{p) 
and dimension Sn{p) 



\/ (^n(p)) = e^^*"^"" 



.Sn{p)/2 



&n(p) 



{R») 



Sn(p)/2 



Sn (p) = 2c(/)p Sz ~ 16y'c(f>p 



(4) 
(3) where the first equality defines the configurational entropy 
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Fig. 4: << -Rg >>g collapses onto a single curve for various 
pressures. Again we show data for our largest system size of 
900 grains. The black solid line is the best fit of a quadratic 
polynomial to the collapsed data. 



of the stress would take an exponentially long time near 
unjamming [B]. We are exploring possible ways of mea- 
suring the mosaic length scale. 
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of force networks, S{p), which can be shown to increase lin- 
early for large p and not to have any extrema or inflection 
points. The comparison in Fig. [3] is based on V {5n{p)) 

obtained from Eq. 01 The quantity B?„{p) 






used in the definition of entropy to remove the force scale 
from the volume of the solution space and assure that the 
entropy is nonncgative. We extract E?{pq) directly from 
the numerics, and have found that it is proportional to 

Conclusions. — We have demonstrated that a static 
correlation function exhibits a diverging length scale as 
the pressure goes to zero in frictionless granular pack- 
ings. This length scale is associated with the PTS correla- 
tion function rather than any two-point correlation func- 
tion. The form of the PTS correlation function reflects 
the structure of the solution space of force networks, and 
can be obtained from the volume of this solution space. 
The constraint satisfaction problem for frictionless disks 
is similar for frictional disks, where now there are addi- 
tional constraints due to torque balance, and we expect to 
see a length scale in frictional materials also. 

In supercooled liquids, the PTS correlation function has 
been associated with the mosaic length scale [21]. In the 
mosaic picture, the supercooled liquid is made up of a 
mosaic of different metastable states. One could specu- 
late that in the granular systems approaching the unjam- 
ming transition, there is such a mosaic characterized by 
the length scale po consisting of regions that are locally 
isostatic, but with different force network solutions. The 
mosaic picture would suggest that the amorphous geome- 
try is made up of distinct subregions, and that stress re- 
laxation would involve rearrangements of the mosaic tiles 
that overcome an entropic barrier which scale as po to 
some power. The time scale for stress relaxation should, 
therefore, increase exponentially as the pressure decreases. 
In the context of a sheared granular system, for example, 
where the shearing is driving the rearrangment of the dis- 
tinct subregions, one would expect that the equilibration 
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